Controlled dynamics of qubits in the presence of decoherence 
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An exactly solvable model for the decoherence of one and two-qubit states interacting with a spin- 
bath, in the presence of a time-dependent magnetic field is studied. The magnetic field is static along 
z direction and oscillatory in the transverse plane. The transition probability and Rabi oscillations 
between the spin-states of a single qubit is shown to depend on the size of bath, the distribution of 
qubit-bath couplings and the initial bath polarization. In contrast to the fast Gaussian decay for 
short times, the polarization of the qubit shows an oscillatory power-law decay for long times. The 
loss of entanglement for the maximally entangled two-qubit states, can be controlled by tuning the 
frequency of the rotating field. The decay rates of entanglement and purity for all the Bell-states 
are same when the qubits are non-interacting, and different when they are interacting. 

PACS numbers: 



I. INTRODUCTION 

Solid-state spin systems have emerged to be promis- 
ing candidates for quantum computation [H, 0, H, 3 and 
communication 0, Q , as they offer the possibility for tun- 
able device parameters and scalability. However, their 
experimental realization requires control over the strong 
environmental decoherence. In most of the solid-state 
systems the interaction of qubits with the environmental 
spins is recognized to be a major source of decoherence 
[7|, l8j . For example in the case of quantum dots the hy- 
perfine interaction of the electron-spin with the nuclear 
spins is the dominant mechanism of decoherence [H, Q . 
In the case of SQUIDS the macroscopic quantum coher- 
ence is shown to be strongly suppressed because by a 
spin-bath consisting of nuclear paramagnetic spins and 
charge defects 10]. 

A well-studied model for decoherence is a central-spin 
system coupled to N non- interacting spins [{j, [Til. [l2l fl3l. 
LLJ) LLy, [lg, Ell ■ An exact evaluation of dynamics for the 
central-spins is obtained for special cases, where there is 
uniform coupling with the bath-spins [IJ, [H, [H|, or a 
special choice of initial states 18], or the system-bath in- 
teraction is a simple interaction between the z component 
of the spins [lj]. Each such cases studied could explain 
the experimentally observed results in some solid-state 
devices [§, [r| . In most of these studies the central-spin 
system either experiences a large static field or zero field. 
But, an important requirement for quantum computing 
is control, which can be either optic al fl9l |2C| or through 
the time-dependent magnetic fields[2l|. In the case of 
single qubits, the control field is an oscillatory magnetic 
field, which can preform gate operations involving arbi- 
trary spin rotations. On the other hand for two-qubits 
the exchange interaction between the qubits is the con- 
trol field, for the implementation of the controlled not 
gate. Recently, coherent spin rotations of a single elec- 



tron were demonstrated in a double quantum dot[2l|. 
Petta et.al., [|[ have implemented the CNOT gate or the 
\J SWAP gate, where the exchange interaction is con- 
trolled by gate voltages. In addition to the control of de- 
coherence of the quantum state, it is equally important 
to control the loss of entanglement in multi-qubit sys- 
tems. For example a maximally entangled state shared 
between two remote sites looses its entanglement due to 
the local environments 16]. In such cases, it is impor- 
tant to know whether a local control could slow-down 
the entanglement loss, thereby increasing the fidelity for 
teleportation. 

In this work we shall address the problem of decoher- 
ence and entanglement loss in the presence of a control 
field: a time-dependent magnetic field. We solve the 
dynamics exactly for the cases of a central-spin system 
consisting of one and two qubits. In the case of single- 
qubit the focus will be on the decay of Rabi-oscillations, 
whereas for the two-qubits it would be on the loss of en- 
tanglement. 



II. MODEL 

We consider a central spin system consisting of one or 
more interacting qubits, coupled to a bath (environment) 
of N spin-1/2 particles. The Hamiltonian for this system 
is given by 

H = [^ S i + (Sfe™* + $~e-*"*)] 



(i) 
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The first two terms represent the self Hamiltonian of the 
system spins {Si), and the last term describes the system- 
bath interaction, where the bath-spins are labeled by /j. 
The qubits are coupled to a common external magnetic 
field B = [B\ cos tut, B\ sin u>t, Bq], and wo,i = eBo,i/ m i 
represents the strengths of the fields along the longitudi- 
nal and transverse directions respectively. The transverse 
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field (rotating field) oscillates with a frequency u>. The 
exchange interaction between the qubits is represented 
by Jij and the qubit-bath couplings by gij . We have ne- 
glected the self Hamiltonian for the bath, as it is assumed 
to be weak in comparison to that of the other interactions 
in the model. Such an assumption is valid, for example, 
in the case of quantum dots, where the qubits are electron 
spins and the bath particles are spin-1/2 nuclei. Since 
the magnetic moment of electron, is much large O(10 3 ) in 
comparison to that of the nucleon, the interaction among 
the bath-spins and their coupling to external fields have 
a comparatively small effect on the dynamics. It was ar- 
gued earlier by Koppens et.al., |22|, that for sufficiently 
large external magnetic fields i.e., uo ^> ^Zj9ij^i, the 
system-bath interaction considered in the above Hamilto- 
nian can replace the usual hyperfine interaction between 
the electron and nuclear spins. Hence the above Hamil- 
tonian can describe the dynamics of the electronic spin 
states in quantum dots, in the regime of strong magnetic 
fields. 

The purpose of the work is two fold, to study the ef- 
fect of a time-dependent field on the decoherence and 
the entanglement loss of the qubits. First in Sec-II, we 
consider the case of a single central spin, and study the 
effect of bath-interaction on the Rabi oscillations between 
the states | f) and | J.). In detail we shall consider two- 
different distributions for the coupling strengths, gu, viz., 
(i) gk — g, V k and (ii) g^ = gexp(— ah 2 ). The Gaussian 
distribution for coupling strengths is of relevance in the 
context of a quantum dot [23(. We also consider the ef- 
fect of bath-polarizations on the qubit dynamics. In Sec- 
III we study the cases of interacting and non-interacting 
spins coupled to the bath. An additional feature that 
the two-spin states have is the entanglement. We shall 
address the question of whether the loss of entanglement 
between the system spins can be controlled the oscillat- 
ing frequency. For both the cases, results are obtained 
analytically. 



III. SINGLE SPIN DYNAMICS AND RABI 
OSCILLATIONS 

In the case of a single central spin the Hamiltonian 
given in Eq[T]gets simplified to 

H = huj a S z + Hljx (S + e~ luJt + S~e luJt ) + S Z ^ 9klk- 

k 

(2) 

The initial state of the qubit-bath system is uncorrelated 
and given by p(0) = ps{0) ® Pb{0), where ps,PB rep- 
resent the density matrices of the qubit and the bath 
respectively. Throughout this paper we shall consider 
the initial state of the qubit to be ps(0) = | T)(T I an d 
the bath state to be 



Pb(0) = PB l ®Pb 2 
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PBi = -x + p B ifm 



The bath spins are polarized along the z direction, and 
the value of their common polarization is denoted by Pb ■ 
For Pb = ±1, p B is pure, and completely unpolarized 
[p B = £ri), for P B = 0. 

The evolution of the total-system is governed by the 
dynamical equation 



P (t) = U(t)p s (0)®p B (0)U\t), 



(4) 



where U{t) = cxp (itTi/h) is the time-evolution operator. 
The reduced matrix of the qubit spin is obtained by trac- 
ing over the bath degrees of freedom, ps{t) = Tr B p(t). 
Before a general expression of U(t) for the full Hamilto- 
nian is given, we shall consider few simple cases. 



A. Zero coupling, g k = 

This is the simplest and well known case studied in 
the context of Rabi oscillations. Setting the system-bath 
interaction to zero, we are left with only the first term 
in the Hamiltonian (see EqJ5]). The unitary operator is 
straightforward to calculate and is given by 



U q (t) = V 



cos ^Q, X 



2i- 



■ {A S Z + mS*} 



(5) 



In the above V is the transformation matrix connecting 
the lab frame and the rotating frame, given by 



V 



e iut/2 Q 
Q-iijJt/2 



(6) 



The Rabi frequency ft — y 'lo\ + Aq, and A = u — ujq. 
The Rabi oscillations occur when Ao = 0, and the tran- 
sition probability for the qubit state to be | \) (J, | at any 
later time is given by 

n(t, W )^ai PS (t)i4>-^|sin 2 ^. (?) 

The frequency of Rabi oscillations depend on the strength 
of the transverse field (coi). 



B. One spin bath 

This is the minimum dimension of the bath, and the 
simplest case of qubit-bath interaction. In the Hamilto- 
nian (see Eqj2]) g^ = 0,V k ^= 1. The unitary operator for 
this case is also straightforward to calculate and is given 
by 



U(t) 



uUt)\ T)b(T 



u-(t)\ l) b (i 



(8) 



where 



U±(t) = V 



t sin -|f2± 

cos -n± X + 2i -| {A±5 2 
2 il± 



.(9) 
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FIG. 1: (Color online) Probability distribution of the resonant 
frequencies given in Eq |13l The bath consists of N = 20 spins. 
The distributions for the qubit-bath couplings are normalized 
such that ^2 k gk = g, where g/h = 20MHz. In plotting the 
above we have set lv = u)a. 



The coefficients Ah 



CO 



(lj q ± gi/2), and Q ± 



A 2 ± . 



Depending on which state the bath spin is, the reso- 
nance can occur either at to = + or lo = ujQ — gi/2. 
For example when the initial state of the bath spin is 
\i/j)b = -^[| T) + I I)], the total state of the qubit-bath 
system can be written as 



\i>) = a + (t)\ H) + b+(t)\ |t) + o_(t)| T4) + b-(t)\ ||). 



(10) 



The cofScients 

a± 
b± 



Jut/2 



cos -f2± 
2 



A± . t 
i—— sin — 11+ 

n 2 ± 



ic 



9. 



(11) 



The above state can be entangled if |a+6_ — 0-b+| ^ 0. 
Since the time-dependent coefficients depend on lo, their 
entanglement can be controlled by tuning the oscillat- 
ing frequency lo appropriately. For u> — loq — gi/2 and 
gi 3> u>i, then the above state at t = n/toi, is close to 
the maximally entangled state = -^=[\ |t ) + | Tl)] as 
a+ = 0, and 6-^0 (ui <C fi-). On the other hand 
when u) = lj + gi/2, the qubit-bath state is close to the 
maximally entangled state \ip) = i[| ft) + | II)]- Note 
that the two maximally entangled states can be differ- 
entiated, by the knowing the value of lj. Even though 
this could be used as a controlled entangling operator, 
or for the Bell-state discrimination, the resonance condi- 
tions are not good from the view point of the qubit, as it 
looses all its polarization. 



FIG. 2: (Color online) Transition probability as a function of 
time. The initial and final states for the transition are | j) 
and | |). The bath is unpolarized (Pg = 0) consisting of 
N — 20 spins. Two different distributions for the qubit-bath 
couplings are considered. The Gaussian distribution for gt is 
normalized such that J^ fe gt = g- In plotting the above we 
have taken ujo = 10 3 , u)\ = 10 and g/h — 40, which are in the 
units of MHz. 



C. N spin bath 

We now go to the more general case of the bath. Since 
the Hamiltonian commutes with the total I z value of the 
bath, its is block diagonal in the I z basis. Hence, all 
it requires is to diagonalize a 2 x 2 matrix of the qubit 
Hamiltonian corresponding to each I z state. The unitary 
operator can then be simply written as 



u{t) = Y J u i q m\ 



(12) 



where \i) are the basis states of the bath and the expres- 
sion for U l q is same as that given in Eq. [9l except that the 
± sign is replaced by the index i. The Rabi frequency 
Qi = ^Jlo\ 
where 



A?, 



depends explicitly on the bath state 



A,- 



LO 



U 



k 



gkik 



(13) 



The frequency at which the resonance occurs is no more a 
single value but rather a distribution determined by the 
coupling strengths g^. For example when all = .9, 
the frequencies Aj are distributed between —gN/2 < 
A, < gN/2. In the limit of large N, this is a Gaus- 
sian distribution centered around A, = lj — ljq. Hence, 
the dominant frequency at which the resonance occurs 
is still located at to — loq. In Fig.l we have plotted the 
distribution of frequencies ^(Aj) for two different dis- 
tributions of the coupling strengths, (i)gk = <7,V k, and 
(ii)(?fc ~ g exp(— 0.01/c 2 ). In both the cases, the distribu- 
tion is centered around uiq, implying that the resonance 
condition for Rabi oscillations is still lo — loq. 
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FIG. 3: (Color online) The 2-component of qubit polarization 
with time. The plot shows the decay of Rabi oscillations with 
time. The bath is unpolarized consisting of N = 20 spins. 
Two different distributions for the qubit-bath couplings are 
considered, where gu — g for both the distributions. In 
plotting the above we have taken wo = 10 3 , wi = 10 and 
g/h — 40, which are in the units of MHz 



FIG. 4: (Color online) Qubit polarization with time. The 
initial bath state is unpolarized consisting of N = 20 spins. 
Two different distributions for the qubit-bath couplings are 
considered, where J^ fc gu — g for both the distributions. In 
plotting the above we have taken uo ~ 10 3 , uji = 10 and 
g/h — 40, which are in the units of mHz. The decoherence 
of one-qubit states (1 — Ps(t) 2 ), can be inferred directly from 
the above figure. 



The reduced density matrix of the qubit can be evalu- 
ated from EqfJ] by tracing over the bath degrees of free- 
dom. Using the above obtained unitary operator, it can 
be written as 



Ps (t) =J2u iPs (0)u}{i\p B (p)\i 



(14) 



For the initial bath state given in Eq|3J the above form 
gets simplified to 



(15) 



where A 4 = (1 + p)*/2+™i(l - p)N/2- m% and m% = 
2(i\ X^fc-^feK)- For the initial state of the qubit ps(0) = 
| T)(T |, the transition probability to the spin down state 
is given by 

p&u) a \ PS (t)\ t) = ^ E A *i§ sin2 ^ ( 16 ) 



In Fig. 2 we have plotted the transition probability 
P±(u>), fixing the time at which the first maximum occurs 
i.e., t = tt/loi. The bath is completely unpolarized con- 
sisting of TV — 20 spins. Two different distributions for 
the coupling strengths gt are considered. Eventhough the 
peak of Pi(u>) occurs at uj — u)q, the width of the distri- 
bution Pj.(w) depends upon the number of bath spins and 
the distribution of g' k s. Either with the increase of the 
bath-size increases or the qubit-bath coupling strength, 
the width of Fj increases, and the probability for the 
transition to the down state decreases. Similarly if the 



distribution g^ becomes narrower, the transition prob- 
ability decreases, and reaches a minimum of 1/2, when 
only one value of g^ is non-zero. To see the decay of 
Rabi-oscillations we evaluate the qubit-polarizations as a 
function of time fixing the r.f frequency at w =wo- 

Using the expression for the reduced density matrix of 
the qubit (see Eq |2"8)) , it is straight forward to evaluate 
the time-dependent polarizations Pg{t) = 2Tr(ps(t)S l ), 
given by 



m) = -^i>[i-2/?(t)] 



Ps(t) 
P V s(t) 



i=l 
2 N 



—— cos uot sm — i h + sin ujt cos — i L 

a 2 2 



cos uot cos — ih — — — sin sin — ih 
2 Qi 2 



(17) 



where fi(t) — ^sin|ilj. We now consider the cases of 
polarized and unpolarized baths separately. 



1. Unpolarized bath 

If we set Pb = in EqJJH then the resultant bath 
state is completely unpolarized, i.e., pu = -^wi- Since 
the state has no polarizations Aj = 1. In Fig. 3 we have 
plotted the variation of z component of the polarization 
with time, for the r.f frequency uj, set equal to ojq. The 
polarization has a damping behavior with time. For the 
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case of uniform coupling strengths gk — g, the above 
expression for the z component of the qubit-polarization 
gets simplified to 



pm = 



i 

2N 



N/2 

E 

=-N/2 



(18) 



For large N, C% /2 _J2 N » exp(-2m 2 /iV) (central 

limit theorem) . Replacing the summation by integration 
over the m values the expression for P§(t) given above 
reduces to 



pm 



cos (Liit + | tan 1 juitj 



(1 



+7 1 



. 7 2 t 2)l/4 

cos (y}\t + I tan -1 r yu>t) 

(l + 7 2 i 2 ) 3 / 4 



(19) 



where 7 = Ng 2 /Aui 2 . One can see that there is an oscil- 
latory power-law decay for long times. The behavior for 
P z (t) would be similar even when the distribution of 's 
are from a Gaussian distribution. The polarization Pg (t) 
plotted in Fig. 4, gives the direct measure of decoherence 
1 - P s {t) 2 for the qubit. 



2. Polarized bath 

We now switch on the bath polarization. For Pb = 
±1, the bath state is pure and the frequency at which 
the resonance occurs get shifted by ±gN/2. For other 
values of polarization (\Pb \ < 1) the bath state is mixed 
and it is difficult to estimate the exact value of resonant 
frequency for a arbitrary distribution of g' k s. In Fig. 5 we 
have plotted the transition probability as a function of 
ui, when the coupling between the qubit and bath-spins 
are same i.e., gk = g- The shift in the resonant frequency 
goes as gPN/2. This shift can be understood from the 
distribution of TOj values. For non-zero Pb, and in the 
large N limit, the eigenvalues of the bath operator I z are 
distributed as 



P(m) = 



TTiV(l-Pl) 



2(m- NPb/2) 2 



Pi) 



(20) 



As the bath polarization increases, the oscillations of 
Pg{t) becomes more stabilized between ±1, indicating 
perfect Rabi oscillations. But the resonant frequency is 
shifted away from ujq by gNPs/2, 



IV. TWO-SPIN DYNAMICS AND 
CONTROLLED ENTANGLEMENT 

In the section we shall consider the case of a central 
spin system consisting of two spin- 1/2 particles which 
can be either interacting or non-interacting. We consider 
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FIG. 5: (Color online) The transition probability as a function 
of r.f frequency to. Initial bath states with different polariza- 
tions are considered. The resonant frequency varies with the 
bath polarization as u) = ujo + gN Pb /2. In plotting the above 
we have taken TV = 20, uio = 100, Ui — 10 and g/h = 1. The 
frequencies are in the units of MHz. 



both the cases. The Hamiltonian describing the two- 
qubit dynamics is given by (see EqfJJ, 



sr 2 e— ')] 



+ JS 1 ■ S 2 + VG^S? + 9 rsi)i z k . (21) 



(2), 



fc 



where S12 = Si + S2 represents the total spin of qubits. 
Depending upon the common set of bath spins with 
which the qubits are interacting, we can either have a 
common or separate baths for the qubits. The total 
spin of the two-qubit system S12 = Si + S2 is invari- 
ant under time-evolution, only in the case of common 
spin environment. Hence for separate baths, we have 
transitions between the singlet and triplet sectors, which 
makes the dynamics more interesting. In case of quan- 
tum dots the overlap between the spatial wave functions 
of the electrons, determines the common and separate nu- 
clear environments for the electrons. For strong overlap 
indicating a large value of J, both the electrons interact 
with the same set of neighboring nuclear spins, and on 
the other hand, when the two electrons are physically 
apart, the overlap is quite small (J small) indicating sep- 
arate nuclear environments. Since J can be controlled 
by gate voltages one can tune from one regime to the 
other. In the case of two-qubits the main interest is 
to study the, the swap operation between qubit states 
(I T-L) I IT))) an d the loss of entanglement for max- 
imally entangled states. Because of the presence of a 
time-dependent transverse field the swap operation can 
controlled both by tuning the rotating frequency u> and 
the exchange interaction J. We use the concurrence mea- 
sure C(t) [2J|, to evaluate the entanglement of the two- 
qubit state. Similar to the considerations made in Sec-I, 
we take the initial qubit-bath state to a direct product 
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p(0) = pi 2 (0) ® Pb(0), where p\ 2 can in general, be rep- 
resented as 



P12 



m,n— 1 



(22) 



where the vector polarizations are given by Pi j2 = 
2Tr[pi2/Si ! 2], and the cartesian components of the tensor 
polarization are U mn = ATr\pi 2 S"F S%]. We note that 
for a pure state (p 2 2 = P12), we have Pi = P 2 < 1, and 
Pi 2 + P t 2 + E m „( nmn ) 2 = 3. The state of the bath is 
same as given in EqJ3J 

We now evaluate the reduced density matrix of the 
two-qubit state, for the cases of common and separate 
nuclear baths. 



A. Common spin-bath for the qubits 

In this section we shall consider the case of a com- 
mon spin environment for the qubits i.e., — g^ in 
the Eqj21j Since, the total spin of the two-qubit system 
is conserved, the triplet and singlet sectors evolve sepa- 
rately. Similar to the analysis for single-qubit we start 
with no bath coupling and then switch on the qubit-bath 
interaction to study the reduced dynamics of the qubits. 
This would help to clearly identify the additional features 
that the bath-interaction can bring to the time-evolution 
of two qubit states. 



1. Zero bath coupling 

When there is no coupling with the bath, the time- 
evolution operator can be found by diagonalizing the 
Hamiltonian in the four dimensional space, by making a 
transformation to the rotating frame. The time-evolution 
operator can then be constructed straight forwardly giv- 
ing 



U l2 {t) = V 12 e 



Jt/i 



+cosAi-j> 2 )(j/ 2 | + \u 3 )(v 3 \} 



+ sinAt{|i/ 2 )(^ 2 |- hX^all + e-^KX^I 



.(23) 



where A = a/Aq + ojf and Ao = u> — u>o- The eigen 
vectors, \i>i) and the transformation matrix are given in 
Appendix A. The above form gets simplified at resonance 
i.e., for lu = ujq. The unitary operator can then be writ- 
ten, in terms of the two-qubit spin operators as 



U 12 (t) = V 12 e lJt ' A 



( -1 + Si ■ S 2 - 2Sf£%) 



+ cosu 1 t(-X + 2SfS%) + isiniv 1 t{Sf + Sf) 



+e~ ut ^l-S l - S 2 ) 



(24) 



Using the above the transition probabilites between the 
states I |t) *->■ | ||) and | fl) <-* | IT) can be easily 
evaluated. They are given by 



P u (t) = ^(l-cos^i) 2 . 



1 , 
2I 1 



AJt 



COS Wit 



(25) 



The transitions within the triplet sector i.e., between 
the states | tt) — ► I II) 1 are controlled only by the fre- 
quency of the transverse field lo\. On the otherhand 
for I IJ,) — ► I |t), which invovles both the singlet and 
the triplet states, the transitions occur both due to the 
exchange interaction and the time-dependent field. For 
u>i =/= J there are two frequencies for Rabi oscillations and 
this would lead to a beat pattern for Pjj (t) . Note that the 
swapping between the qubits in the case of | ||) — * | IT): 
can be controlled both by J and to, which can be useful 
during gate operations. 

It is easy to show that in the transition from | TT) ~~ * 
I ||), no entanglement is generated with time. On, the 
contrary for the transition, | ||) - > | IT)) entanglement 
is generated and it is given by 



C(t) 



1 



1 



e lJt cos 2 u>it 



sin u>it\ 



(26) 



In obtaining the above we have considered the resonant 
frequency for oj. The entanglement here is mainly due to 
the exchange interaction between the qubits i.e, if J = 0, 
C(t) = 0. The time-dependent field, can only tune its 
value, such that the time periods, for which the state 
becomes maximally entangled and unentangled can be 
controlled. We now solve the dynamics in the presence 
of bath interaction. 



2. Coupling to N-spin bath 

As shown in Sec-II that the Hamiltonian has a block 
diagonal structure in the I z basis of the bath, the unitary 
operator here can be written using the above obtained 
form of U\ 2 (t) (see EqJ231 as 



U(t) = YU{ 2 \i)(i\, 



(27) 



where |i) are the basis states of the bath and the expres- 
sion for U± 2 is same as that given in EqJMl except that 
the Rabi frequency A gets replaced by A 



where A; 



uj - 



Efc 9kP k 



A 2 , 



The reduced density matrix for the two-qubit system 
can be written as 



where = 



PS® = wE a ^m°)^ 

i=l 

(1 + P) JV / 2 + m i(l _ p^N/2-m 



(28) 



where m, 
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FIG. 6: (Color online) Concurrence with time. The time- 
dependent concurrence, for the Bell-states \Bi) = 4=[| || 

) + I It)] \B 2 ) = ^[| TT> + I U)],\B 3 ) = ^[| tt) - I II 

)], is plotted for two values of Aw = lu — cjo- The loss of 
entanglement is slowest at the resonant frequency 8u> = 0, 
for all the states, and it increases as the r.f frequency shifts 
away from Wo- The bath is unpolarized with uniform coupling 
strengths between the qubits and bath-spins. In plotting the 
above we have taken N = 20, lo = 100, w\ = 10 and g/h = 1. 
The frequencies are in the units of MHz. 



FIG. 7: (Color online) Concurrence with time. The time- 
dependent concurrence for an initial singlet shared between 
two non-interacting qubits is plotted for three different cases, 
(i) when the r.f frequencies are tuned to their resonance values 
i.e., <5o/ 1,2 ' = 0, (ii) the r.f frequencies are slightly away from 
resonance and (iii) much away from their resonant frequencies. 
In the above Jo/ 1 ' 2 ) = cJq 1 ' 2 ^ — u/ 1 ' 2 ', where uJq 1,2 \ are the 
fields at the local sites of the qubits. The initial states of 
the baths are unpolarized each consiting of N — 20 spins 
and the local fields experienced by the qubits are tjg 1 ' = 100, 



X 1 ) 



= 110 respectively. 



In the presence of the bath interaction, the width of the 
distribution for transition probabilities considered ear- 
lier in EqJ^S] increases, and the amplitude of the Rabi- 
oscillations between the two-qubit states decreases. Even 
these features are same as found in the case of a sin- 
gle qubit, the interest here is the entanglement and its 
variation with the r.f frequency u> and the bath polar- 
ization Pg. It is important to note that the bath does 
not induce any entanglement between the initially un- 
cntangled states, but if there is any non-zero entangle- 
ment at t = it decays with time. For the two qubits 
evolving under their self Hamiltonian (g' k s = 0), there 
is no decay of entanglement for any of the maximally 
entangled states (Bell-states) \Bi) = 4=[| U) ~ \ IT)] 

\B 2 ) = ^[| Tl> + | it>],|Sa) = ^[| TT) + I II)], even- 
though they evolve with time. Because of the bath in- 
teraction the two-qubit state becomes mixed with time, 
and the initial entanglement is no more preserved. 

In Fig. 6 we have plotted the variation of concurrence 
with time for the entangled states mentioned above. The 
rate of loss of entanglement is less, when the r.f frequency 
is tuned to its resonance value (too), and the decay in- 
creases as it is tuned away from resonance. The dynam- 
ics is able to differentiate the evolution of all the Bell- 
states. For non-zero bath polarization the resonant fre- 
quency shifts away from loq by gPsN/2, and this in turn 
controls the entanglement loss. For example if g = 1, 
P B = 1/2 and N = 20, the plots of C{t) for a given 
Bell-state in Fig. 6 would have been reversed as the con- 
dition for resonant frequency has now become 5u> = 10. 



Even though, not shown here, the concurrence and the 
decoherence have similar initial decay rates. Thereafter 
the concurrence falls off rapidly to zero. Similar trends 
are reported (l6j even with a more general Hamiltonian, 
for zero external fields. Note that neither the concur- 
rence nor the decoherence have any dependence on the 
exchange interaction J, as all the states considered above 
belong to a single spin-sector (triplet). 



B. Qubits interacting with separate spin-baths 

In this section we consider the other extreme, where 
the bath-spins with which the qubits are interacting, are 
different. This corresponds to a very weak overlap be- 
tween the electronic wave-functions in the case of a dou- 
ble dot system. Here, the qubit-bath interaction term 
does not commute with the exchange term. This causes 
transitions between the singlet and triplet subspaces, due 
to which the singlet state also decohers. The value of 
exchange interaction will be small in comparison to the 
common spin-bath case, as J is a direct measure of the 
overlap function. We first consider the evolution entan- 
gled states shared between two remote (non-interacting) 
qubits and then switch on J to see its effect on the evo- 
lution. 
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1. J=0 

When there is no exchange interaction, the two qubits 
evolve independent of each other, and an initial direct 
product state will evolve into a direct product state for 
later times, i.e. if ps(0) = Pi(0) ® ps(0), then we have 
ps (t) = pi (t) £g> P2 (t) . Hence the dynamics is of interest 
only when the initial state is entangled. This can be a 
situation where an entangled state is shared between two 
remote qubits and each of the qubit state is decohering 
because of their interaction with their individual baths. 

For a Bell state we have Pi = P% = 0, (see Eql22|) and 
only the diagonal components of II to be non vanishing. 
The two-qubit reduces to form 



P12 



(29) 



Since the time-evolution operator for the present case is 
just the direct product of that of individual qubits, the 
time-dependent polarizations 

2 N 

ir™(t) = * Y^\^\pimsmsmi m 



i=i 



where 



sr(t)sm = E^(*)^E ? ?S ) (*) 5 2- (si) 



The time-dependent matrices yy^ 1 ' 2 ) are given in Ap- 
pendix B. Note that all components of the tensor polar- 
ization become non-zero at a later time. In the absence 
of external fields, only the diagonal components change 
in time, and hence the state at any later time can still 
be written in the above given form (see Eq |2"9"|) . Due 
to this, the expressions for concurrence and decoherence 
at a later time, can be easily obtained in terms of the 
tensor polarizations 16]. But, in the present case, since 
all H mn (t) are non-zero, an expression for concurrence 
C(t) in terms of the tensor polarizations is not straight 
forward. 

In Fig. 7 we have plotted the concurrence when the two- 
qubit state is initially a singlet. We have considered dif- 
ferent cases, where the r.f frequency at each qubit site 
has different values. One can see that for the resonant 
frequencies at local sites the loss of entanglement is min- 
imum. As the r.f frequency moves away from resonance 
the entanglement falls off rapidly. Eventhough, the ini- 
tial loss is same for all the cases, the time at which the 
initial singlet completely looses it entanglement depends 
on the local oscillating frequency. For frequencies much 
away from resonance, there is a limiting curve, as shown 
in Fig. 7. All other Bell-states show the same behavior in- 
dicating that there is no special choice for the maximally 
entangled states when the two-qubits are non-interacting. 

For intermediate regimes where the exchange interac- 
tion is neither strong nor weak, we have few number of 



bath-spins which are common to both and the baths can 
be still be assumed to be separate. As the triplet-singlet 
separation is decided by the value of J, the effect of non- 
zero J for separate baths is to protect the singlet state 
from the effects of decoherence. Any non-zero J would 
differentiate the behavior of C(t) for all the Bell-states. 



V. CONCLUSION 

In this work we have studied an exactly solvable model 
for the dynamics of one and two-qubit states interacting 
with a spin-bath under the influence of a time-dependent 
magnetic field. The magnetic field is static along z direc- 
tion and is oscillatory (r.f) in the transverse plane (x—y). 
The expressions for the reduced density matrices of the 
qubit states at a later time have been obtained. Var- 
ious distributions for the qubit-bath couplings are con- 
sidered. In the case of single-qubit the transition prob- 
ability and Rabi-oscillations between any two orthogo- 
nal states, depends on the strength and distribution of 
qubit-bath couplings and initial bath polarization. The 
strength of qubit-bath interaction, determines the peak 
value and width of the distribution for transition proba- 
bility, whereas the resonant frequency for the transition 
is decided by the bath-polarization. There is a slower 
power-law decay (i 3 / 2 ) for the Rabi-oscillations for long 
times, in contrast to a fast Gaussian decay for small 
times. In the case of two interacting qubits, the pres- 
ence of a r.f field, gives an extra tunable parameter in 
addition to the exchange interaction, for the swapping of 
quantum states between the qubits. Under collective de- 
coherence (common spin-bath) the loss of entanglement 
is different for all the Bell-states, whereas, for the two 
qubits which are non-interacting and decohering due to 
their individual environments, all Bell-states decoher at 
the same rate. By tuning the r.f frequency the rate of 
loss of entanglement can be controlled in both the cases. 



APPENDIX A: EIGEN STATES OF THE 
TWO-QUBIT HAMILTONIAN 

The eigen vectors for the two-qubit Hamiltonian in the 
absence of bath interaction (see Eq |22"]) are given by 

h) = |i)t-— |o) t -|-i)t 

Ml 

M = .I^A^W-^^IO^-I-I). 



A + A s 

\v 3 ) = | I - 2A — ] |1) 7 

l"4> = |0) S 



Ml 

V2(A + A) 



|0) T -|-1) T 
(Al) 



where A = -\/A 2 + lo\ and A = lu — to Q . The above 
eigenvectors are unnormalized and written in the total- 
spin basis, where the basis vectors |1)t = | TT), | — 1)t = 
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II), |0h 



^[| U) + UT)L and |0) s = ^ 
)]. The above eigen- vectors are used in constructing the 
unitary operator Eql23l 

The transformation matrix V12 (see Eq l23[) connecting 
the rotating frame and lab frame is given by 



[| TJ.) - I IT qubits given above. 



V 12 = 



i tuJt 

10 

e~ iwt 

1 



(A2) 



The matrix is written in the total spin basis of the two- 



APPENDIX B: TIME-DEPENDENT MATRICES 

(v (1 ' 2) ) 



The time-dependent coefRcents J7^](i) appearing in 
Eg J31I are given by 



v n? 



- A 2 cos Qit) — Ajfij sin flrf — 2A, 



Aitti sin Qit 
2A i u> 1 sin 2 |Qj 



rif cos Hit 

LOi^li sin f2,t 



Wi sin 2 |f2j 
sin f2jf 
(Af +lo 2 coaflA) 



(Bl) 



In a similar 

,( 2 ) 



where fij = -\/A 2 + w 2 and A. 

way the matrix ?)( 2 ) can be written with Aj = uj^ - . 
Here ' ; represents the local fields at the qubit sites 
and a/ 1 ' 2 ), the respective local r.f frequencies. The trans- 
formation matrix R = cos(u)^i)X+ 2i sin(a;W . 
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